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The Vaisala-Brunt frequency is examined for data obtained at an 
anchor station in the Florida Current for a ten day period. There 
are ninety-three separate STD stations examined to yield a mean 
profile of the Vaisa la-Brunt frequency at this location. 

The maximum frequency of the mean profile is calculated to be 
0.0169 sec-1 and opcurs at a mean depth of 100 meters. This cor- 
responds to a Vaisala period of 6.19 minutes. 

The depth at which the maximum frequency occurs varies with 
time^ but with no apparent periodicity; this is also true of the 
magnitude of the maximum frequency. This randomness may be the 
predominant feature of the Vaisala-Brunt frequency. 

The plot of frequency versus depth for each station reveals fine 
structure in the vertical column^ showing a layered structure. The 
shallower layers appear to have a stronger lamination effect on the 
density structure. 
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1. INTRODUCTION 

Density gradients affect the motion of the seas* The salinity, 
temperature and pressure gradients all cause the sea water density to 
vary with depth* Most of the density gradient is caused by the compres- 
sion of the water^ the rest by the variation of temperature and salinity 
with depth (Eckart, 1961). Associated with the density gradients are 
well defined density boundaries which exhibit certain frequencies and 
modes of oscillation. Thus, any given density boundary will have its 
own normal oscillating frequency* This is called the Vaisala-Brunt 
frequency and is expressed in radians per second. It is a measure of 
the stability of oceanic stratification (LaFond and LaFond, 1966). 

The distribution of this frequency is one of the most important dynam- 
ical features of the ocean. 

The concept of vertical stability in the ocean was expressed by 

Sverdrup, ^ in 1942 as E = 10’^ dcr t , and appears in Richardson's 

dZ 

number as 1 d [/^ ] . Pollack (1954) defines vertical stability as 
[yo ] dZ 

follows : 

1. As a measure of the restoring force (or resistance to change) 
in the vertical direction. 

2. As an acceleration per unit change of geometrical depth; the 
acceleration may be expressed per unit change of pressure 
also; Hesselberg used this approach. 

3. As an expression of any continuous function which is a 
parameter of the resistance to deformation of the fluid column. 





^s. 
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4. As the amount of energy required to move a particle of unit 
mass a unit displacement. 

In definitions 2 and 3^ stability is a function of the difference 

between an observed vertical mass distribution and that of a neutral 

equilibrium mass distribution. Therefore^ zero corresponds to the 

neutral position^ stability is indicated by the positive values and 

instability by negative values (Poliak^ 1954). 

According to Hesselberg^ stability £ = g ^ y or E = g [ d/^ - (d/^ ]y 

r dZ r [dZ (dZ)^] 

where (d/) is the density gradient of neutral equilibrium or the mass 
(dZ)^ 

distribution in a column of homogeneous water (pressure=0^ salinity and 
temperature=cons tants ) with negligible molecular conductivity. Bjerknes^ 
^ (1934) noticed the adiabatic density gradient { dP ) was the re- 

(dP )a 

ciprocal of the term in the Lapalcian expression for the velocity of 
2 

sound C = (dP) ^ (Poliak, 1954). Therefore, the adiabatic density 

2 

gradient can be replaced by 1/C (see Appendix A). This is possible 

because of the analogy between the pressure, p, and the geometric unit 

of depth, Z. Poliak equates stability in the geometric units as: 

E= 9.81 [ dP - 9. 81 P ] , which appears similar to an expression devised 
P [dP ~c^ ] 

by the meteorologist Vaisala in 1928, which he designated N^, a sta- 
bility expression with dimensions of frequency. Eckart mentioned this 
expression in notes of a seminar given in 1949. According to Poliak, 
the source for the expression was a lecture by Solberg at the Univer- 
sity of California in 1946. The original reference could not be located. 

We can examine the frequency, N, by noticing that if the density 
of a fluid is changed adiabatically, then dP = Cp dP . VTtere dP and dP 
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are the changes in pressure and density^ and C is the speed of sound. 
If dP is caused by changing the level of a parcel of water by an 
amount dZ then dP = -^gdZ if we assume a positive Z - axis up^ and g 
is the acceleration of gravity. 

Elimination of the term dP gives: 

i ^ 

dZ ^ 

The quantity g/C^ is identified as the adiabatic density gradient and 
is that part of the actual density gradient associated with compres- 
sibility effects. Let E = - (1 d/^ ^ find that^ generally 

(7® dZ c2) 

speaking, it is a positive value in the seas. The quantity N in- 
troduced by Vaisala and Brunt is equal to gE. It has frequency 

dimensions and therefore the period, T = ^ (Eckart, 1961). 

N 

In considering the significance of the coefficient N, Eckart 

(1960) considers the oscillation of a small fluid mass displaced from 

its zero position and allowed free, adiabatic movement, with the 

positive Z-axis upward. He assumes the particle of fluid is in a 

small balloon- like membrane which permits changes in density to occur 

unimpeded, and thereby maintain the outside and inside pressures the 

same. The displacement causes a pressure change and an accompanying 

density change within the balloon, so that d/o inside = The out- 

C2 

side density change will be : 

dP outside = ^ where r is the displacement. The buoyant 
^ dZ ^ 

force on this balloon is then equal to 



g (d /^outside - d inside), (1) 
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which is 

g ^ (df - ^g) =-P n2 i, (2) 

^ (dZ 

where 

= -g [i ^ + S ] (3) 

[/° dZ C2] ’ 

If the frequency^ is real^ the parcel of water in the balloon would 

oscillate freely at the frequency N. Instability would exist if 
was negative^ causing the ballon to leave its initial position (Eckart, 
1960). 

It is believed that N is the limiting frequency for the free 
internal waves because of viscous damping (McLellan^ 1965). Internal 
waves can exist if their frequency of oscillation is greater than the 
inertial frequency^ but less than the maximum Vaisala frequency^ 

^max‘ There should^ therefore, be no free internal waves with fre- 
quencies above N , and few with frequencies above the local Vaisala 
frequency in a particular water mass (Cox, 1962). In examining in- 
ternal waves near Bermuda, Haurwitz, ^ £l* ; (1959) used the Vaisala 
frequency. They associated the maximum frequency with the seasonal 
thermocline and the second maximum found at a deeper depth, with the 
permanent thermocline. They also concluded that at any depth, there 
should be no frequency larger than the maximum Vaisala frequency. 

Hence, this may be designated the ’’cut-off” frequency. This cut-off 
frequency is known to vary with time and space and must be computed 
for each vertical density distribution. When converted to period the 
Vaisala frequency relates to the natural oscillation governed by the 
stability of the water layers (Lee, 1960). 

In the upper four thousand meters of the ocean the period has a 
wide range, from one minute to four hours. In deeper waters the 
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frequency is difficult to compute because the terms in the brackets 
of the equation are very nearly equal and of opposite sign. At the 
temperatures and pressures found in the deep waters the values of the 
coefficient of thermal expansion are doubtful, and at best are extra- 
polated values. These extrapolated values create large uncertainties 
in the value of Vaisala frequency. Uncertainty also arises from 
errors in measuring small salinity gradients (Eckart, 1961). 

An investigation of the Vaisala frequency in the Florida Current 
was conducted using pilot data collected in a twenty-four hour period 
during the winter of 1966 by Dr. Saul Broida- Two locations were used 
with stations taken alternately between them. Density parameters 
were measured U\ s 1 1 u with a Bisset-Berman Salinity-Temperature-Depth 
(STD) system. The data were analyzed at fifteen meter depth increments 
from the surface to seven hundred meters. The Vaisala frequency and 
period were computed for each depth increment^ and plotted versus 
depth. 

The frequency equation as stated by Eckart was used: 

CL d/^ + g ) in radians per second. 
dZ C2) 

A computer program was devised to calculate the terms of the equation 

and solve for the frequency and period, where T = 2Tr . Curve fit- 

N 

ting was accomplished by a cubic spline Interpolation (see Appendix B). 

Analysis of the resulting pilot data disclosed some interesting 
features : 

1. At any given depth the frequency varied with time. Some 
periodicity was implied, but insufficient data points did 
not allow any definite conclusions. 
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1 . Minimum period values occurred at the boundary of the 
isothermal layer and the mixed upper layer. 

3. The shape of the frequency curves indicated a multi-layered 
structure of laminations of density, as noticed elsewhere 
in the oceans by Stommel and Fedorov (1967), and others. 

4. The depth of maximum frequency varied with time. 
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II. METHODS 



A. DATA COLLECTION 

During t)ie period from 19 February to March I, 1968 the R/V 
PILLSBURY was anchored in the Florida Current in 460 fathoms of water 
at Latitude 25-39. 1 N Longitude 079-48.0 W (Fig. 1). The objective 
was to gather temperature, salinity, and depth information to study 
density fluctuations in the Florida Current, and analyze the Valsala- 
Brunt frequency at this location. This cruise, P6801, was the longest 
continuous anchor station completed in the Florida Straits. 

A Bissett-Berman Sal Inlty-Temperature-Depth probe (STD) system 
was utilized to gather the in situ data. Probe casts were conducted 
at two hour intervals. The initial lowerings were made to 750 meters, 
but an unfortunate wire entanglement later limited the casts to a 
maximum 500 meter depth. The STD probe utilizes sensing elements which 
send back the Jji situ data as a frequency signal for each parameter. 
These frequencies are recorded on graphic plots and magnetic tape. 

The conversion to actual values of depth, temperature, and salinity 
is done with standard conversion tables of frequency versus the para- 
meter measured. To maintain the calibration and accuracy of the 
system, Niskln bottle samples were taken for comparison. At selected 
times a sample was taken at the surface and one at the maximum depth 
along witlr the probe lowering. The salinity and temperature readings of 
both Instruments were compared and found to be in good agreement. 
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Because of the nature of the Florida Current a strong surface 
current prevailed throughout the experiment, approximately 4,2 knots. 
Technical problems and rough seas limited the total number of stations 
taken to ninety- three during the ten day period. 

The ship's position was determined by visual bearings and radar 
ranges to Fowey Rocks Light. The accuracy of the anchor position is 
+ 100 meters in an East-West direction; + 400 meters North-South. 

Some gaps do exist in the sequence of data collection which are 
attributed to mechanical failures and weather conditions. However, 
the major portion of the stations are at two hour intervals with the 
last sequence of stations, number seventy-one to ninety- three, being 
taken at hourly intervals. 
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FIGURE 1. 



Chart of ship's position 
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B. DATA ANALYSIS 

The STD data from the cruise were in analog format for each 
station, but those stations taken after station number 32 also recorded 
the parameters on magnetic tape via a digital data log unit (DDL). 

Those stations not recorded on magnetic tape initially were digitized 
with a CALMA converter. This is a manual trace of the temperature 
and salinity graphs, which records the horizontal and vertical trace 
distances on magnetic tape. By computer conversion the distances are 
changed to the desired parameters « The end result is a magnetic tape 
record of all the station data. 

The STD scanning rate was set at 0.1 seconds giving a high density 
of points on the traces and tapes. To select the most representative 
value of a parameter at a given depth an averaging process was used. 

The DDL values were averaged at twelve meter intervals from the sur- 
face to the maximum depth of the station cast. A comparison between 
the averaged values and the actual graphic values of temperature and 
salinity at given depths indicated a tolerance of + 0.02^^ C for tempe- 
rature and + 0.03 o/oo for salinity. The averaging process provided 
the best values for analysis for the total cast. 

The Temperature-Salinity (T-S) diagram for the composite of the 
data was investigated and compared to the historic data for that region. 
The historic data was taken from the works of Parr and various cruises 
conducted by the Institute of Marine Sciences. Even though the curves 
were geometrically coincident, the salinity values of some stations 
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appeared to be lower than the historic data. Although the spread of 
the T-S curve was larger than might be expected^ the relative changes 
in temperature and salinity appear adequate for analysis of the gra- 
dients. Thus, these relative values are acceptable for the investi- 
gation of the Vaisala frequency which is density gradient dependent. 

There are many small features shown on the temperature traces 
which are probably real, but the salinity traces have many small scale 
spikes of low salinity which appear false, probably caused by imperfect 
matching of the time constants of the temperature and conductivity 
sensors of the STD. The falsification appears as fine scale variations 
in salinity at the points of strongest vertical temperature gradient. 
The deeper the traces go, the more the salinity trace becomes a sharp 
indicator of temperature gradients. 

It is felt that many of these spikes are compensated for in the 
averaging process and also in the interpolation of points while curve 



fitting. 
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C. COMPUTATION 

Once the station data was compiled in terms of temperature, 
salinity, sigma - t, and specific volume anomaly, for each twelve meter 
increment of depth, the parameters for solving the Vaisala-Brunt fre- 
quency equation were avialable. The value for the acceleration of 
gravity was assumed constant at 9.7893 meters per second at Latitude 
25.7 N. Sound speed was computed at each depth increment by the 
following equation: 

C = 141,000 + 421 t - 3.7 + 110 S + 0.018 d, in (cm/sec), (4) 

where 

t = temperature (°C), S = salinity (o/oo), d = depth (cm). 

These values of speed were converted to meters per second. 

The i^ situ density was computed and used in determining the 

quantity dP . The equation used to find the density was: 
dZ 

OC (35,0,P)+S = OC situ (cm'Vqm), (5) 

where 

OC(35,0,P)= standard specific volume at depth, and S = specific 
volume anomaly. The reciprocal of the specific volume ^ situ is the 
in s itu density, . The standard specific volume was computed by 
linear interpolation. The values decreased linearly with depth by 
45 X 10"^ cm^/ gm/ meter. 

The values of density situ were plotted versus depth for each 
station and a cubic spline interpolation fit was made to the points. 
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(Fig. 2). Using interpolated values from these curves the slopes, 

dP , were found at five meter intervals, as were the values of density 
dZ 

at those depths. 

The frequency equation was solved using these interpolated values 
in the equation: 

= -8 (1 ( 6 ) 

(r dZ c2) 

Frequency was converted to period by T = 2 TT minutes. The frequency 

N(60) 

was plotted versus depth and analyzed, (Fig. 3). A mean frequency 
curve was also plotted for the entire cruise data by averaging the 
interpolated values of frequency at five meter increments for the 
ninety- three stations (Fig. 4). Where negative values of N were 
encountered, they were eliminated from the average to present only 
positive values of frequency. These negative values are interesting 
in that if they really exist they could represent instability in the 
laminar structure of the water column, but they may be caused by 
erroneous data readings as a result of time lags in the sensors. It 
would be interesting to investigate these in future studies. The 
historic data from the Florida Current indicate strong stability in 
the area of the experiment; therefore, it is felt that these negative 
values might really be neutral stability points, approaching N=0, 

Because they were few and variable in depth, it is felt that the 
averaging process was not prejudiced by their elimination, nor was the 
investigation of the frequency in depth and time. V/hen possible, depths 
were chosen where a minimum or no negative values of N were noted 
during analysis. 

To help diminish errors in using Vaisala frequency as a measure 
of stability, small depth increments were chosen for analyzing gradients 
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FIGURE 2. 



Density versus depth profiles 



(6 stations) 
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FIGURE 3. 



Frequency versus depth profiles 



(6 stations) 



I 






18 





u 







i 



19 



FIGURE 4. 



Mean frequency versus depth profile 
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in density. Hence^ the five meter increment for interpolated values 
of density. According to Roden (1968) it is believed that any fun- 
damental oscillations would differ from the Vaisala frequency by about 
twenty-five per cent for this vertical displacements In studying the 
Vaisala frequency in the Gulf of Alaska^ Roden found that the non- 
linear cases showed stability oscillations no longer symmetrical about 
the mean, as is the Vaisala oscillation, and the fundamental fre- 
quency of oscillation differed from the Vaisala frequency. 

In computing the standard deviation and variance of the fre- 
quency for the entire period of observations, 100, 245, and 350 meters 
were selected as representative depths. They were so chosen because 

o 

they contain a minimum of negative values, which were to be ex- 
cluded; they included the mean depth of the maximum frequency as 
compiled in the mean frequency curve; and they are spaced far enough 
apart to give an in depth view of the frequency on a statistical 
basis . 

From each station the maximum frequency was gathered and plotted 
versus time (Fig. 5), and so was the depth at which the maximum fre- 
quency occurred (Fig. 6). These were visually inspected for some 
recognizable periodicity, but none was readily apparent. In fact, 
the random nature of the function is the primary observation. 

From the mean profile of Vaisala frequency a maximum value of 
0.0169 sec"^ was obtained, which corresponds to a period of 6,19 
minutes. The depth of the maximum frequency was 100 meters. The 
computed standard deviation of depth for maximum frequency was jh 15.2 
meters, with a relative dispersion (coefficient of variance) of 16.6 



per cent. 
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The standard deviation of the frequency at 100 meters was + 
0.0037 sec’^^ with a relative dispersion of 22.8 per cent. At 245 
meters the standard deviation was + 0.0030 sec‘^ and the relative 
dispersion was 40.5 per cent. For 350 meters the standard deviation 
was + 0.0031 sec"^ with a relative dispersion of 60.9 per cent- 
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FIGURE 5. 



Maximum frequency versus time. 
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FIGURE 6. 



Depth of maximum frequency 



versus time 
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III. RESULTS 



A. DISCUSSION 

Results of the experiment discussed in this thesis lead to the 
observation that the Vaisala-Brunt frequency in the Florida Current 
near the axis of the current is a complex feature of stability with 
no apparent periodicity with respect to depth or time- An examination 
of each plot of frequency versus depth shows that the maximum fre- 
quency occurs in the thermocline where the mean density gradient is 
greatest and decreases above and below this level as the water becomes 
more homogeneous. The smallest frequency values are observed in the 
deeper layers and near the surface. 

The weather during the data collection phase was particularly 
stormy at the anchor station. This may account for the well mixed 
upper layer and fairly abrupt thermoclines recorded. Thus^ we ex- 
perienced a three layer oceanic structure^ with relatively homo- 
geneous water above and below the thermocline layer. The plots of 
frequency versus depth show a sharp maximum frequency at or just below 
the thermocline boundary with smaller values of N elsewhere. The 
thermocline thickness appears to be between 100-150 meters. 

It is observed that the STD data provide continuous profiles of 
the thermohaline^ and therefore density^ structure. The instrument 
makes obsolete the older methods of interpolating between widely spaced 
sample points. The STD reports discontinuities varying in magnitude^ 
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time^ and space as the rule rather than the exception. In the plots 
of frequency versus depth (Fig. 3)^ series of layers are observed 
clearly with sharp boundaries separating each from its neighboring 
layers, and these exist throughout the water column. The gradients 
or boundaries between the layers are often thin and sharp. The 
shallower layers show stronger lamination effects in the density 
structure than do the deeper depths. These are obvious observations 
of the data collected and are only the beginning of exploration into 
the recognition of this fine structure. Stommel and Fedorov (1967) 
present some interesting ideas on the dynamic processes involved in 
this fine structure and possible gradient effects on vertical mixing 
and mean turbulent transport. Data similar to that obtained for this 
thesis could be analyzed in the light of their studies and add much 
to the study of the Florida Current. 

It is recognized that the Vaisala-Brunt frequency plays a major 
role in the study of internal waves, therefore its distribution with 
depth in the Florida Current is of interest. As noted^ it is a very 
complicated function of depth. It appears that the complex structure 
occurs because of temperature and salinity perturbat ions . 

The average plot of frequency versus depth presents a representa- 
tive trace of the Vaisala-Brunt frequency for the Florida Current 
during the period of the cruise. The small deviations indicate that 
the average frequency figure might be useful in future calculations. 

B. RECOMMENDATIONS 

In looking ahead to future work, it is recommended that this type 
cruise be conducted again in the summer months, preferably to obtain 
a seasonal variation of the density structure and see the effect on 
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the Vaisala-Brunt frequency. More emphasis should be placed into 
examining the fine structure, however, because this may Lead to a 
better understanding of how the density structure varies, and why<, 

It may also offer an understanding of the Vaisala frequency randomness 
of variation. 

Project MIMl, a study of acoustic signal propagation across the 
Florida Straits, collected acoustic data s imulaneous ly with our col- 
lection of density data. An examination of these data might provide 
an insight into the internal wave configurations of the Florida Current, 
or produce a means of measuring Vaisala-Brunt frequency directly, if 
a definite relationship exists between the variations of N and the 
acoustic variations received in Bimini. It offers an excellent op- 
portunity to compare physical parameters with acoustic data taken 
simultaneously over an extended period, along the acoustic axis. 

Any future investigations of the Vaisala-Brunt frequency in the 
Florida Current and its distribution will enhance the dynamical 
investigations of this region. 
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IV. SUMMARY 

The Vaisala-Brunt frequency was calculated and investigated for 
a given location in the Florida Current near the axis of the current, 

The analysis of data from ninety-three separate STD stations gives a 
composite profile of the Vaisala-Brunt frequency versus depth for this 
region- The experiment was designed to obtain density information for 
an investigation of the fluctuations of the frequency as a function 
of time and depth. The results of the investigation are station 
profiles of the frequency versus depths a mean profile for the aggregate 
of stations^ a plot of maximum frequency versus time and the depth of 
N max versus time. The standard deviation and variance were computed 
at three selected depths for frequency^ and for the depth of N max. 

The investigation reveals that the frequency is a very complex 
function of depth and time which appears to vary in time^ space^ and 
magnitude in a random manner; typical of many oceanic processes. This 
randomness of the function is the most significant observation, A 
representative figure for N max is presented along with the depth, at 
which it occurs for possible use in future work related to internal 
waves^ acoustics, or density stability in the Florida Current, 

CONCLUSIONS 

The mean value of maximum Vaisala-Brunt frequency at this location 
in the Florida Current is 0,0169 sec’^ with a corresponding period of 



6 . 19 minutes . 
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The average depth of occurrence of the maximum Vaisala-Brunt 
frequency is 100 meters with a standard deviation of + 15.2 meters 
and a coefficient of variance of 16.6 per cent^ or relative dispersion. 

The standard deviation of the mean Vaisala-Brunt frequency at 
100 meters is + 0.0037 sec’^ with a coefficient of variance of 22.8 
per cent. At 245 meters the standard deviation decreases to + 0.0030 
sec*^ and the coefficient of variance increases to 40.5 per cent. At 
350 meters the standard deviation is + 0.0031 sec*^ with a coefficient 
of variance of 60.9 per cent. All the coefficients of variance are 
increasing with depth, with the standard deviation fairly constant with 
depth. The deviation is slightly higher within the thermocline, near 
maximum frequency. 

The average plot of frequency versus depth is shown in figure 4 
and is representative of the Vaisala frequency in depth for this 
location. 

The depth at which the maximum frequency occurs varies with time 
as does the magnitude of maximum frequency, however, no apparent 
periodicity is noticed. The randomness of this feature may be its 
most important characteristic. 
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APPENDIX A. VELOCITY OF SOUND TERM IN THE 



VAISALA- BRUNT EQUATION 
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From Hook's Law: STRESS = K X STRAIN. 

E = Bulk Modulus, in dimensions of pressure FORCE or MASS x ACCELERATION 



area 



AREA 



/° = density, in dimensions of MASS 

VOLUME 



so that 



E = ML 

Wl.2 



M 

lt2 



and /° = M , 



1? 



E=M xL^ = l2=(L)2 

/° LT? M 'P' ( T ) 



( 1 ) 

( 2 ) 



which is in velocity units squared, 



According to Newton, the speed of propagation of waves 

C^ = E , 

r 



( 3 ) 



substituting pressure P, for E, we obtain the same results, so that 



P = 
7 



(^) 



3 



and for an adiabatic change 

(dP ) = C2. (5) 

)a 
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APPENDIX B. 



CUBIC SPLINE INTERPOLATION 



1 
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A new approach based on an old idea is used to interpolate a set 
of data. It is based on the idea of using a different expression in 
each sub- interval ^i+l^ auxiliary conditions of smoothness 

at the joints instead of a simple closed form expression^ valid over 
the entire interval. It was developed because the other methods (e.g,^ 
least squares) failed to give good results under certain conditions; 
at the two extremes of a set of data or could not provide a good 
evaluation of the first derivative along the computed polynomial fit. 

The basis for all interpolation schemes is that given a discrete 
set of n + 1 points an nth degree polynomial can be passed. An inter- 
polating polynomial of degree as high as the number of data points 
would allow can be obtained. But this can be impractical sometimes. 
Some of the problems encountered with interpolating schemes like the 
least square method are: 

!• A set of data having points which are far apart from 
each other^ will fail to give a good approximating 
polynomial, since the conditions for smoothness are 
very fragile, in most cases. 

2. If the change in the slopes from point to point are 
very sharp, the resultant polynomial produces gross 
errors, since this rate of change of slopes has a 
direct effect, not only where the change occurs but 
also along the entire interval. 

3. Theory states that the higher the degree of computed 
polynomial the better the approximation. Again, 
experience shows that one could go only so high on the 
degree of the interpolated polynomial before the errors 
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pile up making the whole thing frustrating, when a 
considerable amount of data points are available. 

Experience shows that polynomials higher than 13th 
degree are impractical. But this depends on the 
distribution of the data points. 

All this led to the development of the Spline interpola t ion which 
is based on the idea of passing a low degree polynomial through a few 
points at a time making use of conditions of continuity and existence 
of first and second derivative- 

A cubic polynomial is interpolated in each sub interval [X|“l,Xj^], 
i = l,...n. It is accomplished with the following conditions met, 

1. Continuity of the function 

2. Continuity of first derivative 

3. Continuity of second derivative. 

Due to the nature of the computat ions , data which is equally spaced 

provides short cuts to the procedure. With this in mind the data for 
this experiment was chosen at equal intervals of depth. 

For a mathematical treatment of the Cubic Spline Interpolation 
scheme the reader is referred to *'Approximat ion of Functions” by 
Carl de Boor (1965). Elsevier Publishing Co., New York. 
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APPENDIX C. 



STATISTICAL FORMULAS 
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mean frequency X = ^ (1) 







n 










n = 


number of points 




STANDARD DEVIATION 


s = 


jd - 


(X)^ 


(2) 




\ 


n 






VARIANCE 


V = 






(3) 


COEFFICIENT OF VARIANCE 


(relative dis 


per Sion) 






C = 


s = 


standard deviation 


(4) 






X 


mean 
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